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We propose angle-resolved photoelectron spectroscopy of aerosol particles as an alternative way 
to determine the electron mean free path of low energy electrons in solid and liquid materials. The 
mean free path is obtained from fits of simulated photoemission images to experimental ones over 
a broad range of different aerosol particle sizes. The principal advantage of the aerosol approach 
is twofold. Firstly, aerosol photoemission studies can be performed for many different materials, 
including liquids. Secondly, the size-dependent anisotropy of the photoelectrons can be exploited 
in addition to size-dependent changes in their kinetic energy. These finite size effects depend in 
different ways on the mean free path and thus provide more information on the mean free path than 
corresponding liquid jet, thin film, or bulk data. The present contribution is a proof of principle 
employing a simple model for the photoemission of electrons and preliminary experimental data for 


potassium chloride aerosol particles. 

PACS numbers: 34.80.-i, 79.60.-i, 82.80.Pv, 82.70.Rr 

I. INTRODUCTION 

The inelastic and elastic mean free path (MFP) of an 
electron is the average distance an electron of a given ki¬ 
netic energy travels in a material between two successive 
inelastic or elastic collisions, respectively [1, 2]. Electron 
mean free paths (EMFPs) are important quantities for 
a number of physico-chemical phenomena. These range 
from aerosol physics and chemistry, which influence the 
radiative energy balance of the Earth and other planets, 
to photoemission processes in interstellar dust clouds, ra¬ 
diotherapy, and nuclear plant operation (see [3-18] and 
references therein). A complete understanding of the 
whole photoemission process from the absorption of light 
to the emission of electrons from the material is crucial 
in this context [19-21], For many compounds, the values 
of the inelastic and elastic MFP are not known or in¬ 
accurate, largely due to various experimental difficulties. 
Especially challenging are experiments at electron kinetic 
energies (KEs) below 100 eV and for liquids, which can¬ 
not be brought into vacuum as bulk [22-26]. While sev¬ 
eral methods have been introduced for the energy range 
from 15 eV to 100 eV, the uncertainties of the mean free 
paths remain high in the range below 15 eV [26, 27]. 
Recently, progress has also been made for liquids by em¬ 
ploying liquid microjets [28, 29]. 

Here, we introduce an alternative way to determine 
the MFPs of slow electrons with KEs below 20 eV which 
relies on angle-resolved photoelectron spectroscopy of 
aerosol particles. Aerosol particles are ionized by single 
photons from a vacuum ultraviolet (VUV) light source. 
The three-dimensional (3D) distribution of the gener¬ 
ated photoelectrons is then projected with a velocity map 
imaging (VMI) setup [30] on a two-dimensional (2D) 
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position sensitive electron detector [10, 31] (section III). 
To characterize the properties of photoelectrons, we ex¬ 
tract four parameters from the photoelectron images in 
addition to the electron kinetic energy: the total elec¬ 
tron yield and three anisotropy parameters, which de¬ 
scribe the anisotropy of the photoelectrons with respect 
to the propagation and the polarization direction of the 
ionizing radiation. These four parameters depend in a 
characteristic way on the EMFPs and the aerosol parti¬ 
cle size. The EMFPs are extracted from fits of the sim¬ 
ulated anisotropy parameters to the experimental ones 
varying the EMFPs. The simulated photoemission im¬ 
ages are obtained from a simple model (section II). Sim¬ 
ilar models have been widely used in the past, also for 
particulate matter (see [32-38] and references therein). 
Briefly, classical electrodynamics is employed to treat 
light-particle interaction and a simple three-step model 
is used to describe electron formation, transport in, and 
emission from the aerosol particles. 

The use of the aerosol approach has two main advan¬ 
tages compared with other methods. Firstly, it allows 
one to determine EMFPs also for liquids, which can¬ 
not be brought into vacuum as thin films or bulk due 
to their high vapor pressure. This is not a fundamen¬ 
tal problem for aerosol droplets. Secondly, aerosol data 
recorded over a range of different aerosol particle sizes in 
principle contain more information on the EMFPs than 
corresponding thin film, liquid jet, or bulk data, because 
the size-dependence of the anisotropy can be exploited 
in addition to the kinetic energy. The key point is as 
follows: The variation of the aerosol particle size changes 
the total path length of the ionizing radiation and the 
total scattering path lengths of the electrons inside the 
particles. In addition, the variation of size corresponds 
to a variation in curvature. This changes the geometry 
of electron escape from the particle. More importantly, 
the curvature characteristically modifies the pattern of 
the light intensity inside the particle and thus the ioniza- 
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tion probability. In other words, the curvature strongly 
influences the anisotropy of the photoelectrons. Since all 
these finite size effects depend in different ways on the 
EMFP they provide characteristic and complementary 
information on the EMFP. 

The purpose of this contribution is a proof of prin¬ 
ciple of the idea. We use potassium chloride (KC1) as 
a model substance because many of its properties are 
well known and aerosol formation is straightforward [39- 
42]. Furthermore, calculations for EMFPs at low KEs 
already exist for KC1 [43], which allow for a comparison 
with results from the present work. In section IV, we 
first present the main prediction from modeling, which 
we then compare with preliminary experimental data for 
KC1 aerosol particles. 



FIG. 1. (Color online) Calculated intensity patterns of the 
ionizing radiation inside KC1 nanospheres with a radius of 
5 nm (panel a), 100 nm (panel b), and 500 nm (panel c), re¬ 
spectively. The calculations are for a photon energy of 10 eV. 
The light intensity decreases from yellow to red to blue color. 


II. MODEL AND PARAMETER DEFINITION 

The model consists of three parts: The photonic model 
(subsection II A) describes the interaction of the aerosol 
particle with the electromagnetic radiation. The electron 
emission model (subsection IIB) includes the formation 
of the electron at a given location inside the aerosol par¬ 
ticle, the transport of the electron through the particle, 
and its emission at the particle surface. The electron 
projection model (subsection IIC) defines the projection 
of the 3D photoelectron distribution onto a 2D detector 
plane. This subsection also introduces the parameters, 
a, /^backward, /^forward, and 7 , that describe the character¬ 
istics of the photoemission for comparison with experi¬ 
mental data. We use a 3D Cartesian coordinate system 
with the origin at the center of the aerosol particle. The 
photons travel in the positive x direction and are linearly 
polarized along the y axis. The aerosol particles travel in 
the negative y direction and the photoelectrons are pro¬ 
jected in the positive z direction onto the 2D detector 
plane which is parallel to the x-y plane. 

A. Photonic model 


culated with the discrete dipole approximation (DDA) 
using the ADDA package [33, 44-47]. We tested the accu¬ 
racy of this approach with the Lorenz-Mie method which 
provides an analytical solution for spherical particles and 
with previously published data [32, 48-52]. 

Figure 1 shows the squared amplitude of the inter¬ 
nal field \E(a)\ 2 , i.e. the light intensity, in the x-y plane 
through the center of a spherical particle with a radius 
of 5 nm, 100 nm, and 500 nm, respectively, at a photon 
energy of 10 eV. The light intensity inside the smallest 
particles is homogeneous. Only a very minor decrease in 
intensity in the propagation direction of the light (posi¬ 
tive x direction) is visible, which is due to minor absorp¬ 
tions in the particle. The largest particle absorbs the 
light almost completely already in the outermost surface 
layer of the particle that faces the incoming radiation. 
There is no light that penetrates into the particle’s core. 
For intermediate sizes, here a 100 nm particle, the light 
intensity inside the particle is strongly structured with 
regions of low and high intensity. As mentioned, these 
size-dependent intensity patterns together with the in¬ 
fluence of the size-dependent curvature on the electron 
emission determine the anisotropy in the photoelectron 
images. 


Classical continuum electrodynamics is employed for 
modeling light-particle interaction. For the KC1 par¬ 
ticles, we use the frequency-dependent refractive index 
data in the VUV from Ref. [39] and assume spherical 
particles surrounded by vacuum. Note that we find essen¬ 
tially the same results for spherical KC1 particles as for an 
ensemble of randomly oriented KC1 cubes with rounded 
corners, consistent with the results for NaCl in Ref. [38]. 
The ionizing radiation is modeled as a monochromatic 
plane wave which is linearly polarized along the y-axis, 
providing a good model for the synchrotron radiation 
used in the experiment. We assume that the probability 
of generating an electron at a specific location a inside 
the particle is proportional to the intensity of the inter¬ 
nal electric field E(a) at this location (subsection IIB). 
The electric fields inside and outside the particle are cal- 


B. Electron emission model 

We use a simple three-step model for the description 
of the electron emission [35]: The first step is the gener¬ 
ation of the photoelectron at location a inside the par¬ 
ticle. Once formed, it is randomly scattered inside the 
aerosol particle (second step) until it reaches the surface 
and is emitted from the particle into vacuum (third step). 
For the first step, we assume that the probability p(a , 9) 
of generating an electron at location a inside the parti¬ 
cle with momentum in direction 6 is proportional to the 
light intensity (subsection II A) and the classical dipole 
emission probability, i.e. 

p(a,6) (x \E(a)\ 2 cos 2 (9) 


(1) 
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where E{a) is the internal electric field at location a and 9 
is the angle between the directions of the internal electric 
field vector E{a) and the initial velocity vector ?Ti n itiai(d) 
of the electron. The initial electron KEs Ek,\ nitiai are 
chosen randomly from a Gaussian distribution modeled 
after the photoelectron spectrum of gas phase monomer 
KC1 [53] in order not to include scattering a priori. For 
photoelectron spectra recorded at photon energy hv, they 
fulfill the condition 


0 Ek. initial ^ hv Eg 


( 2 ) 


where E g is the band gap, for which we use a value of 
8.4 eV [41]. The Gaussian is centered at hv — E g — A, 
where A = 0.34 eV is the half width at half max¬ 
imum [53]. We have found that the exact shape of 
the curve and the exact initial KE distribution does 
not significantly influence the asymmetry parameters a, 

/^backward ■ and pforward ■ 

The second step, the migration of the electron within 
the particle, is modeled as a random walk with the inelas¬ 
tic MFP and the elastic MFP as step lengths for inelas¬ 
tic and elastic scattering events, respectively. A random 
walk corresponds to isotropic scattering of the photoelec¬ 
trons in the particle. This is probably a good assump¬ 
tion here because the directionality of scattering in bulk 
KC1 is weak [43] and because our KC1 particles are likely 
polycrystalline. The probabilities of inelastic scattering 
events p\ and elastic scattering events p E is calculated as 


Pi = l/(s + 1) and p E = s/(s + 1) . 


(3) 


where s is the ratio of the inelastic and the elastic MFP. 
We use s = 0.1 in this work consistent with Ref. [43]. As 
both elastic and inelastic scattering events are isotropic 
in our model, the value of s does not significantly influ¬ 
ence the anisotropy parameters a , /3backward, and /3forward- 
Their values for a 100 nm radius particle, for example, 
vary by less than 1% if s is varied between 0.1 and 10. 
Note that scattering of low energy electrons in solids, 
in particular elastic scattering, is far from being under¬ 
stood [43]. For the mean energy loss per inelastic scat¬ 
tering event, the energy loss parameter e, we use 

e = 10 meV at 300 K. (4) 


This corresponds to the value 9 meV at 270 K from 
Ref. [43], which has been corrected to account for the 
higher temperature of about 300 K in our experiments. 
For this correction, the temperature dependence of the 
mean phonon energy from Ref. [14] is used. Again, we 
can show that the exact value of e has only a minor in¬ 
fluence on the asymmetry parameters. However, it will 
significantly affect the electron yield. Since we did not 
perform yield measurements in the present work, we use 
10 meV throughout the paper. 

For the third step, the emission of the electron at the 
particle surface into vacuum, the electron must fulfill the 
classical momentum threshold condition for escape 


r < arccos 


Ea/ Ax,surface 


(5) 



FIG. 2. (Color online) Example photoelectron images with 
different anisotropy parameters a, ^backward, and /3f orwar d. 
The panels a and b are representative for the situation found 
for small and large particles, respectively. The backward and 
forward plane is the left and the right half-plane, respectively. 
The electron intensity decreases from yellow to red to blue 
color. 


where r is the angle between the electron velocity vec¬ 
tor and the surface normal and Ek, surface is the elec¬ 
tron KE [37]. Both values are taken just before ejection; 
i.e. at the surface inside the particle. We use a value of 
Ea = 0.5 eV [41] and subtract it from the normal compo¬ 
nent of the electron KE just after ejection into vacuum. 
For every ejected electron, we store its velocity vector 
and its KE resulting in a 3D electron emission pattern. 
The corresponding velocity vector and the corresponding 
KE in the vacuum just after emission are referred to as 
^vacuum and E Kl respectively. 


C. Electron projection, anisotropy parameters, 
photoelectron yield 

Experimentally, it is very challenging to record this 
3D electron emission pattern. Instead, we use velocity 
map imaging (VMI) [30] in our experiments (section III), 
which projects the emitted 3D electron distribution onto 
a 2D position sensitive electron detector with the help of 
an electrostatic lens system. The recorded 2D photoelec¬ 
tron images are referred to as ” experimental raw images” 
in the following. To model the experimental projection, 
we map each electron onto the detector plane with be~ l& , 
where b oc sjE k — Ek, z and 0 is the azimuthal angle in 
the x-y plane with respect to the y axis. The correspond¬ 
ing photoelectron images are referred to as ’’calculated 
raw images”. 

To describe the angular dependence of the electron 
emission patterns, we extract from the raw images three 
anisotropy parameters, denoted as a, /?backward, and 
/^forward- Note that in VMI the 3D emission pattern can 
be retrieved from a reconstruction of the 2D raw image 
only if the electron distribution has an axis of cylindri¬ 
cal symmetry parallel to the detector plane [54]. In our 
case, however, the light attenuation inside larger parti- 
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cles (Fig. 1) breaks this symmetry so that reconstruction 
is not possible. Therefore, we determine the anisotropy 
parameters directly from the calculated or experimen¬ 
tal raw image. For their definition, we split the image 
plane into two half-planes with x < 0 and x > 0 re¬ 
ferred to as ” backward plane” and ” forward plane”, re¬ 
spectively (Fig. 2). The anisotropy parameter a describes 
the anisotropy of the photoelectron image with respect 
to the propagation direction of the light. It is defined as 
in Refs. [38, 55] as 

a = I x+ /I x - ( 6 ) 

where 7 X+ and 7 X _ are the number of photoelectrons pro¬ 
jected onto the x > 0 and x < 0 half-planes, respec¬ 
tively. For gas phase molecules, a = 1 because there is 
no anisotropy in the propagation direction of the light 
for these samples. For KC1 aerosol particles, by contrast, 
one expects a values between 1 and 0 as a result of the 
attenuation of the light by absorption inside the particle 
(Fig. 1). For very small aerosol particles, a is still close 
to 1 because light absorption is minor here (panel a of 
Fig. 2). With increasing particle size, a decreases con¬ 
tinuously. An example for an image with a = 0.25 is 
provided in the panel b of Fig. 2. 

In VMI photoemission imaging, the anisotropy is usu¬ 
ally described by the anisotropy parameter /3, which is 
defined through /(</>) oc [1 + /3/2(3cos 2 </> — 1)] [56-58]. 
The angle (f> is the polar angle with respect to the po¬ 
larization direction of the light (y direction) and I(<j>) is 
the number of photons detected in angle (f) in the recon¬ 
structed image (center slice image). Here, we keep the 
same definition for /?, except that we determine it di¬ 
rectly from the raw image and separately as /^backward in 
the backward half-plane and as /3f or ward in the forward 
half-plane (Eq. 7). As a consequence, these /? param¬ 
eters do not have the same physical meaning as usual. 
Rather they are phenomenological parameters used to 
characterize the anisotropy associated with the polariza¬ 
tion direction of the light. We define them by 

I{(t>~) oc [1 + /?backward/2(3cos 2 </>“ - 1 )] and , , 

I{(t> + ) OC [1 + /3forward/2(3 COS 2 <fr + - 1)] 

where </>“ and </> + are the polar angles in the raw image 
with respect to the polarization direction of the light (y 
direction) in the negative (x < 0 ) and positive (x > 0 ) 
half-planes of the image, respectively, and I(4>~) and 
7(0“) are the numbers of photoelectrons detected in di¬ 
rections (\r and </> + , respectively. Examples of raw im¬ 
ages with different /3 parameters are shown in Fig. 2. 
Although the three asymmetry parameters are not com¬ 
pletely independent, they all contain complementary in¬ 
formation which is useful for the extraction of electron 
MFPs. 

Lastly, we define the total photoemission yield 7 as a 
fourth parameter to characterize the electron emission. 
The parameter 7 is the ratio of the number of photoelec¬ 
trons that escape from the particle g Q ut and the number 


of absorbed photons c7? 2 otaI , i.e. 


7 = 9out/c7? 2 otal = PaE 2 (a)/E 2 otal (8) 


where q ou t = c^2, a p a Ef{a) and c is a measure of the 
absorption cross section. The index a runs over all pho¬ 
toelectrons that are formed at location a in the parti¬ 
cle. Ef(a) is the local field intensity at location a and 
7f 2 ota i = J]] a 7? 2 (a). We assume the probability of gen¬ 
erating an electron at location a to be proportional to 
the local field intensity at this position and p a £ { 0 , 1 }. 
p a is 1 if the electron generated at location a is finally 
emitted from the particle and it is 0 for electrons that 
are generated in the particle but never emitted from the 
particle. 


III. VELOCITY MAP IMAGING EXPERIMENT 

We have performed preliminary photoemission experi¬ 
ments on KC1 aerosol particles with radii of 37.5, 50, 75, 
100, 125, and 150 nm at photon energies of 10, 11.4, and 
12.4 eV at the VUV beamline of the Swiss Light Source 
at the Paul Scherrer Institut. The KC1 aerosol particles 
are generated in an atomizer (TSI Inc., model 3076) and 
dried after the exit of the atomizer with a home-made 
diffusion dryer. After size-selection in a scanning mo¬ 
bility particle sizer (SMPS, TSI Inc., model 3080), they 
enter the photoelectron spectrometer through a home¬ 
made aerodynamic lens [4, 10, 31]. Typical particle den¬ 
sities after size selection vary from 0.6 • 10 s to 3.8 • 10 5 
particles/cnr 3 . VUV synchrotron radiation is used for 
single-photon ionization of the aerosol particles. The 
generated photoelectrons are projected onto a 2D elec¬ 
tron imaging detector (Roentdek, 40 mm diameter) with 
a home-made VMI lens system [10, 31, 59]. Typical ac¬ 
quisition times are on the order of several hours per im¬ 
age. 

High background signal and time constraints at the 
beam line limited the quality of the images and made 
measurement on smaller particles impossible, while fluc¬ 
tuations in the photon flux made photoelectron yields 
unreliable. In the uncertainty estimates given in sec¬ 
tion IV B, we have tried to account for the most sig¬ 
nificant sources of error. Hence, the experimental data 
we present are only preliminary and could be significantly 
improved and extended by further optimization of the ex¬ 
perimental setup by reduction of background noise and 
increased photon flux stability. However, section IV B 
demonstrates that even these preliminary experimental 
data strongly support our proof of principle study. 
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FIG. 3. (Color online) The lines are simulated anisotropy pa¬ 
rameters a, ^backward, and /^forward and electron yield 7 for in¬ 
elastic MPFs of 0.2 (dashed red), 0.5 (solid red), 0.75 (dashed 
green), 1 (solid green), 1.5 (dashed blue), 2.5 (solid blue), 5 
(dashed cyan) and 10 nm (solid cyan). The experimental data 
are shown as black crosses with error bars. All data are for 
KC1 particles photoionized at 10 eV photon energy. Please 
note that the values of /^forward for different MFPs cross each 
other. 


IV. RESULTS AND DISCUSSION 
A. Model predictions 

Figure 3 shows the evolution of the anisotropy param¬ 
eters and the total electron yield as a function of the size 
of the KC1 particles as predicted by the model. The dif¬ 
ferent lines represent values for eight different inelastic 
MFPs from 0.2 to 10 nm (see caption of Fig. 3). For 
these simulations, we assume at first that the inelastic 
MFP is constant during the random walk, i.e. that it does 
not depend on the KE of the electron. Figure 3 reveals 
that all four parameters depend on the inelastic MFP, 
but in different ways. The parameter a depends only 
weakly on the MFP for particle radii below about 25 nm 


because the light intensity in these particles is more or 
less homogeneous (Fig. 1) and the number of scattering 
events is low to moderate. For larger sizes, however, the 
values of a become more clearly distinguishable for dif¬ 
ferent inelastic MFPs, which is mainly a consequence of 
the inhomogeneous light intensity inside the aerosol par¬ 
ticles (Fig. 1). Moreover, 7 shows a similar trend. It is 
less sensitive to the MFP for smaller particles than for 
larger ones, for which it shows a distinct dependence on 
the MFP. Just the opposite behavior is observed for the 
two /3 parameters. In contrast to a and 7 , /^backward and 
/^forward are clearly more sensitive to the MFP for par¬ 
ticle radii below about 50 nm than for larger particles. 
The complementary and characteristic dependence of the 
four parameters on the inelastic MFP is useful in that if 
one covers a broad aerosol particle size range one obtains 
comprehensive information on the MFP. This, in turn, 
may allow one to determine accurate values for MFP us¬ 
ing aerosol measurements. 

For the simulations in Fig. 3, we have assumed that the 
inelastic MFP remains constant throughout the random 
walk of the electron in the particle. In reality, however, 
the electron MFP can pronouncedly depend on the KE of 
the electron, which then changes during the random walk. 
From the observed dependence of the four parameters on 
the MFP in Fig. 3, one would expect to find noticeable 
effects on these parameters when the energy dependence 
of the inelastic MFP is included. Fig. 4 confirms this 
expectation. This figure shows the dependence of the 
four parameters on the inelastic MFP as a function of 
particle size for three different cases: a constant inelastic 
MFP of 1 nm (identical to the solid green line in Fig. 3) 
and two inelastic MFPs that depend in different ways on 
the KE of the electron. The energy dependences of the 
latter are derived from the calculated MFPs for acoustic 
phonon-electron interaction and optical phonon-electron 
interaction provided in Fig. 1 (c) of Ref. [43]. For both 
types of interaction, the energy dependence is to a good 
approximation linear in the region that is of interest here. 
We have approximated them by two straight lines with a 
slope of +0.8 nm/eV and -0.8 nm/eV, respectively, both 
of them sporting the value of 1 nm MFP at 1 eV elec¬ 
tron KE. The influence of the energy dependence of the 
MFPs on the anisotropy parameters in Fig. 4 is in gen¬ 
eral moderate in particular for the smaller particles. An 
exception is /^forward for larger particles, which exhibits 
a pronounced energy dependence in this size range. The 
parameter /3f or ward is mainly determined from electrons 
that originate from the particle half-sphere in which for 
larger particles the intensity pattern of the light varies 
sensitively with particle size (Figs. 1 and 2). As a conse¬ 
quence, /^forward is here also very sensitive to the inelastic 
MFP. The strong dependence of the electron yield 7 on 
the MFP is not a surprise since it scales with the MFP 
(Eq. 8 ). 

The dependence of the four parameters on the refrac¬ 
tive index n + ik is illustrated in Fig. 5, for the exam¬ 
ple of a 50 nm KC1 particle and for a constant inelastic 
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FIG. 4. (Color online) Anisotropy parameters a, /3forward, 
and /^backward and electron yield 7 simulated for a constant 
inelastic MFP of 1 nm (green dots), an energy dependent 
inelastic MFP representing the interaction with the acous¬ 
tic mode (blue triangles), and an energy dependent inelas¬ 
tic MFP representing the interaction with the optical mode 
(red squares) [43]. The experimental data are shown as black 
crosses with error bars. All data are for KC1 particles pho- 
toionized at 10 eV photon energy. 


MFP of 1 nm, again neglecting the dependence of the 
MFP from the electron KE as in Fig. 3. The calcula¬ 
tions are for the KC1 refractive index data between 9 and 
13.5 eV displayed in the two top traces in Fig. 5. The 
green dots represent calculations for which only the re¬ 
fractive index and the wavelengths of the light vary while 
all other parameters are kept constant. In particular, the 
initial KE distribution Ek, initial is not varied according 
to the incident photon energy (section IIB). We chose 
a constant Ek, initial distribution as calculated for a pho¬ 
ton energy of 10 eV (section IIB). The variations in a, 
/^backward, ^forward in Fig. 5 follow more or less the varia¬ 
tions of the refractive index. The changes in the absolute 
values of a, ^backward , and /3f or ward are moderate to pro¬ 
nounced. The absolute value of 7 , by contrast, is hardly 


influenced by the refractive index. A similar behavior is 
also observed for other particle sizes (data not shown). 
For the black squares in Fig. 5, the Ek, initial distribution 
is varied in addition to the refractive index and wave¬ 
length. The Ek, initial distribution at each photon energy 
is calculated as described in section IIB. The comparison 
with the green dots (dashed lines) reveals that the addi¬ 
tional influence of the initial kinetic energy is not very 
pronounced, except for 7 at very low photon energies. 
Fig. 5 reveals that the quality of refractive index data is 
important for the determination of accurate MFPs from 
anisotropy parameters. 


B. Comparison with experiments 

The modeling results clearly reveal the principal ad¬ 
vantage of aerosol studies over other methods for the de¬ 
termination of accurate electron MFPs . The advantage 
lies in the characteristic dependence of the anisotropy pa¬ 
rameters on the EMFP as a function of aerosol particle 
size (Fig. 3). If a broad size range is covered the comple¬ 
mentary information contained in the three anisotropy 
parameters allows one to extract detailed information on 
the EMFP. Given that high quality experimental data are 
available, one should be able to extract accurate EMFPs 
from a fit of calculated to experimental data, includ¬ 
ing information on the energy dependence of the EMFP 
(Figs. 4 and 5). 

To complete our proof of principle study, this section 
provides a comparison of the modeling results with pre¬ 
liminary experimental data for KC1 aerosol particles. As 
mentioned in section III, the quality of the experimental 
data is limited and it was thus not possible to extract reli¬ 
able data for the electron yield 7 or to record meaningful 
data for very small aerosol particles. We have recorded 
photoemission images at photon energies of 10, 11.4, and 
12.4 eV for six different particle sizes. As examples, we 
present mainly the data for 10 eV photon energy and 
only a summary for 11.4 eV photon energy. Because 
of the low signal to noise of the measurements recorded 
at 12.4 eV and overlapping signals from background gas 
phase water we do not provide any data for this photon 
energy. In view of the limited quality of the experimen¬ 
tal data, we do not attempt to extract any energy de¬ 
pendent information on the inelastic MFP (Fig. 4), but 
compare it first only with the model predictions for con¬ 
stant MFPs. This comparison is shown in Fig. 3. Despite 
the preliminary character of the experimental data, both 
the experimental and the calculated anisotropy parame¬ 
ters show the same trends with increasing aerosol particle 
size. The values of a and /^backward decrease continuously 
with increasing particle size while /^forward increases with 
increasing size. All values of a and /3backward are con¬ 
sistent with a short MFP. The comparison of these data 
with the model predictions results in an average inelas¬ 
tic MFP of about 1 nm at an average electron energy of 
about 1 eV. A similar analysis of a and ^backward for the 
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FIG. 5. (Color online) The two top traces show the refractive 
index data n + ik from Ref. [39] as a function of the photon 
energy. The four bottom traces show the anisotropy param¬ 
eters a, /^backward, and /^forward and the total electron yield 7 
as a function of the photon energy. All calculations are for a 
constant inelastic MFP of 1 nm and a KC1 particle radius of 
50 nm. For the green dots, it is assumed that only the refrac¬ 
tive index varies with photon energy. For the black squares 
it is assumed that the refractive index and the initial kinetic 
energy distribution vary with photon energy. Note that the 
lines in the four lower panels are merely meant to guide the 
eye. 


measurement at a photon energy of 11.4 eV results in 
an average inelastic MFP of about 0.5 nm at an average 
electron energy of about 2.5 eV. Both values are in accor¬ 
dance with the theoretical prediction of Ref. [43]. Note 
that our inelastic MFP corresponds to the dominating 
inelastic MFP, i.e. the shorter one in Ref. [43]. On the 
basis of the present data, we are unable to judge how 
significant the agreement with Ref. [43] is. It should be 
noted in this context that Akkerman et al. use a simi¬ 
lar model as the present work; i.e. they also treat the 
electrons as classical particles that move along stochas¬ 
tic trajectories. The main difference to our model lies 


(a) light propagation (b) 



FIG. 6 . (Color online) Experimental (panel a) and calculated 
(panel b) raw photoelectron image of a KC1 particle of 150 nm 
radius. The particles are ionized with light of 10 eV photon 
energy. 


in the MFPs. They use a quantum mechanical model 
to calculate the MFPs while we extract them from a fit 
to experimental data. At least, the comparison suggests 
that the modeling approach chosen allows for a meaning¬ 
ful analysis of the experimental data. 

In the above estimates of the MFPs, we have not used 
/^forward- Fig- 3 shows that /3forward is only consistent 
with an inelastic MFP of about 1 nm for the small par¬ 
ticle size, but not above about 100 nm. The same phe¬ 
nomenon is found for the data for 11.4 eV photon energy 
(not shown). At this point, we can offer two explanations. 
The inconsistency of /3f 0 i- W ard observed for larger particle 
sizes could be a consequence of the low signal to noise in 
the photoelectron images. Figure 6 shows experimental 
and calculated photoelectron images for a 150 nm KC1 
particle. The parameter ^forward is defined in the half¬ 
plane in which the number of detected electrons is low 
for larger particles (Fig. 2) because they originate from 
the particle hemisphere that is only weakly illuminated 
(Fig. 1). Therefore, the signal to noise ratio is lower 
and background contributions in the experimental image 
are higher in the half-plane which is used to determine 
/^forward compared with the one that is used for /3backward- 
This can result in ill-defined values for ^forward- The in¬ 
consistency of /^forward for larger particles might also orig¬ 
inate from the energy dependence of the inelastic MFP. 
As discussed, Fig. 4 shows that /3f or ward is more sensitive 
to this energy dependence compared with ^backward and 
a. The comparison of experimental and calculated data 
in this figure also reveals that the energy dependence of 
the inelastic MFP can easily account for the apparent in¬ 
consistency of /^forward- The data for the interaction with 
the acoustic mode and a 1 nm inelastic MFP do no longer 
show any systematic deviation of /3f or ward for larger parti¬ 
cle sizes. For an energy dependence of the inelastic MFP 
according to the acoustic mode, all three anisotropy pa¬ 
rameters provide consistent results for the inelastic MFP. 
More accurate measurements will resolve which of these 
effects dominates, background contributions or the en¬ 
ergy dependence of the inelastic MFP. But already the 
preliminary experimental data discussed here provide a 














convincing proof of principle for the new aerosol approach 
to determine MFPs of slow electrons. 


V. CONCLUSIONS 

This proof of principle study shows that angle-resolved 
aerosol photoelectron imaging opens up a new and poten¬ 
tially more accurate way to determine electron MFPs for 
low KE electrons in solids and liquids, provided that ac¬ 
curate experimental data and a realistic model are avail¬ 
able. The advantage of the aerosol approach arises from 
the characteristic information on the MFP that is con¬ 
tained in the anisotropy of the photoelectrons. The vari¬ 
ation of the particle size - and thus the variation in cur¬ 
vature - results in characteristic modifications of the light 
absorption and the electron transport inside the aerosol 
particles and thus of the angular dependence of the pho¬ 
toelectrons. Our modeling results clearly demonstrate 
that the use of the anisotropy information in principle 
allows one to extract detailed information on the MFPs 
from fits of simulated photoelectron images to experi¬ 
mental ones. 

Important for the determination of accurate MFPs is 
the quality of the experimental photoelectron images. 
Crucial are a good signal to noise ratio, low background 
contributions, a quality electron imaging system, and 
stable particle and photon fluxes. The first two factors 
strongly influence the quality of the anisotropy data. The 
third factor is crucial to obtain reliable data for the total 
electron yield. For the present study, we use a velocity 
map imaging photoelectron spectrometer. We do not fit 
simulated photoemission images directly to experimental 


ones. Instead, we first define four parameters (the elec¬ 
tron yield and three anisotropy parameters) from the im¬ 
ages, which we then use for the comparison between sim¬ 
ulation and experiment. A major reason for this choice is 
the quality of the preliminary experimental data. In fu¬ 
ture studies, it is planned to directly fit whole images in 
order to exploit their full information content. The elec¬ 
tron MFP is a model-dependent quantity. Here, we used 
a simple conventional model as a start. The incorpora¬ 
tion of more sophisticated models might further increase 
the utility of the aerosol approach [3, 24, 42, 43, 60]. 

The aerosol approach is particularly attractive for very 
low-kinetic energy electrons and for liquids for which ex¬ 
perimental data on MFPs are urgently needed. MFPs of 
low-kinetic energy electrons are important in many dif¬ 
ferent fields ranging from astrophysics to radiotherapy, 
and liquid phase chemistry. We have performed first ex¬ 
periments for various liquids and are currently applying 
the aerosol approach to these systems. 
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